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Abstract 

The Small Subgraph Conditioning Method has been used to study 
the almost sure existence and the asymptotic distribution of the num- 
ber of regular spanning subgraphs of various types in random rf-regular 
graphs. In this paper we use the method to determine the asymptotic 
distribution of the number of 3-star factors in random d-regular graphs 
for d > 4. 
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1 Introduction 

It their remarkable papers (see [TO] and [TT]) Robinson and Wormald showed 
that for d > 3 and dn even, a random d-regular graph contains a Hamilton 
cycle with probability tends to 1 as the number n of vertices tends to infinity. 
They used the Small Subgraph Conditioning Method (see [B] or [H] for details) 
to prove the existence (with high probability) of perfect matching in such 
graphs when n is even. 

The method has been used to determine the existence with high probabili- 
ty of, and the asymptotic distribution of, the number of /c-regular spanning 
subgraphs (for k = 1,2) and the number of long cycles in random (f-regular 
graphs (see [5j, [7j and |9j). 

A star is a tree with at most one vertex whose degree is greater than 1. 
A k-star is a star with k leaves. A k-star factor in a graph is a spanning 
subgraph whose components are /c-stars. 

We use notations P (probability), E (expectation) and Var (variance). We 
say that an event Yn occurs a.a.s (asymptotically almost surely) if 

lim Py„ = 1. 
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In [2] Assiyatun and Wormald have used the method to investigate the 
a.a.s of 3-star factor in random d-regular graphs. This is the first time the 
method apphed to non-regular subgraphs in such graphs. 

Assiyatun and Wormald [2J started proving the existence a.a.s of a 3-star 
factor in random (i-regular graphs for d > 4 by showing the existence a.a.s 
of a 3-star factor in random 4-regular graphs. Then using the contiguity 
of models of random regular graphs (see [12]), they obtained the existence 
a.a.s of a 3-star factor in random (i-regular graphs for fixed d > 4 as desired. 

As a completion to the result in [2j, in this paper we use the method 
to determine the asymptotic distribution of the number of 3-star factors 
in random d-regular graphs for d > 4. However, due to the complexity of 
some part of the computation we are only able to obtain the asymptotic 
distribution for 4 < d < 10. Nevertheless, in most part of the computation 
we obtain the result for general d > 4. The main result obtained in this 
paper is presented in the following theorem. 

Let Qn,d be a probability space contains of d-regular graphs with n ver- 
tices. In asymptotic statements about properties of Gn,d, we restrict n to 
even integers when d is odd. 

Theorem 1.1 Restrict n to mod 4 and 4 < d < 10. Then G G Gn,d cl-cl.s 
has a 3-star factor. Furthermore, letting denote the number of 3-star 
factors in G G Gn,d, 



EYr 



d 



W = J\{1 + Skf'^e'^"^'^ forn^oo 

k=3 

where are independent Poisson variables with EZ^ = \k for k > 3 and 
{d-1)^ 



2k ' 

k 




4{d-l){d-3/2) 



As in [9] we will first work on the pairing model which was first intro- 
duced by Bollobas (see [H]). This model can be described as follows. Let 
V = ljr=i ^ be a fixed set of dn points, where = for every i. A pairing 
is defined as a perfect matching of points of V into dn/2 pairs. A pairing 
P corresponds to a random d-regular pseudograph G(P) in which every Vi 
is regarded as a vertex and each pair is an edge. We use Vn^ to denote the 
probability space of all pairings. 

As shown in |[3], the probability that the pseudograph has no loops or 
multiple edges (i.e simple graph) for a fixed d, is asymptotically bounded 
below by a positive constant. Moreover, each simple graph arises with the 
same probability as G{P) for P £ Vn,d- Hence using the following property 
we obtain the desired result in Gnd- 
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Lemma 1.1 A property of graphs that holds a.a.s for random pseudographs 
arising from Vn,d will also hold a.a.s for Qn,d- 

Given two sequences a„ and hn, we denote a„ ~ 6„ if ^ — 1 for n — )■ oo. 
We denote falling factorial n(n — 1) • • • (n — m + 1) by [ri\m and Stirling's 
formula by 



n! ~ yj2T:n ( — ) for n — )• oo. 



2 The variance of the number of star factors 

Throughout this paper we define 

(2m)! 



N{2m) 



ml T 



as the number of perfect matchings of 2m points. 

Note that counting subgraphs of the pseudograph coming from Vn^d is 
equivalent to counting the corresponding sets of pairs in the pairing. For 
that purpose, parallel edges are distinguishable from each other (especially 
as they come from distinct pairs in the pairing). 

Let n = mod 4 and define as the number of 3-star factors in G{P) 
coming from Vn,d for d > 4. We have the following theorem. 



Theorem 2.1 



d{d-3/2r/^-'/ui.] i '-'ir-^' ] ) . (2.1) 



Proof. See SI for details. □ 



2_y/2 nd-l){d -2)Y^^' 
d'^ [ 3! 



Using the method in P] (see |3j and [U] for similar argument) we obtain 
the following theorem. 

Theorem 2.2 [13\ Restrict n to mod 4 and define as the number of 
3-star factors in G{P) coming from Vn,d for 4 < d < 10. Then 

Proof. We count the ways to lay down an ordered pair of 3-star factors in 
P £ Vn,d. In general, a set of pairs in P inducing a subgraph of a given type 
will be called by the same name in the pairing. 

Let Si be a 3-star factor of P for i = 1,2. Let T = 5i n 5*2. Define 
T = U5=i Tj, where (see Figure 2.1) 

(i.) Ti consists of xi 1-stars, Si and ^2 have only one common leaf. 
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(ii.) T2 consists of x-2 l-stars, and ^2 have one common leaf and one 
common center, 

(iii.) T3 consists of 2-stars, S\ and ^2 have two common leaves, 
(iv.) Ti consists of X4 3-stars, Si S2, 

(v.) consists of X5 0-stars, 5i and ^2 have only one common center. 




Figure 2.1: Two intersecting 3-star factors 
Given Si , the number of possibilities of the intersection T is 



ra/4 

Xl 



^xi - xi^ f "/4 -X1-X2 



X2 



X3 



3x3 



n/4: — Xl — X2 — X3 — X4 



n/4 — Xl — X2 — X3 

X4 

(n/4)! 3^1+^2+3:3 



I X4 



X5 J {n/4: — Xl — X2 — Xs — X4 — x^y. Xi\x2lx3\x4\x5\ 

(2.3) 

There are — xi — X2 — X3 — X4 — x^) 3-stars in both 3-star factors that 
don't share any edge. We call these edge-disjoint 3-stars isolated 3-stars. 

We have to complete 52 by creating the isolated 3-stars and completing 
Ti,T2,T3 and into 3-stars. 

The centers of the isolated 3-stars in S2 can not be chosen from the vertex 
set of T nor the center of the isolated 3-stars in Si. There are (2x1 + 2x2 + 

3X3-|-4x4-|-X5)-|-(n/4-.Ti-.T2-X3-X4-X5) = (n/4 -|- Xi -|- X2 -|- 2x3 -|- 3X4) 

such vertices. Thus the number of ways to choose these centers is 



3n/4 — Xl — X2 — 2x3 — 3x4 
n/4 — Xl — X2 — X3 — X4 — X5 



(3n/4 — Xl — X2 — 2x3 — 3,T4)! 



(n/4 — Xl — X2 — X3 — X4 — X5)! (n/2 — X3 — 2x4 + x^y. 

(2.4) 

The number of ways to choose the leaves of the isolated 3-stars in ^2 is 



n/4— 0:1—2:2—2:3— X4— 2:5 — 1 

n 

fc=0 



/ 3n/4 — Xl — X2 — 2x3 — 3x4 — 3fe 
V 3 



(3n/4 — Xl — X2 — 2x3 — 3x4)! 



(31)71/4-0:1-0:2-0:3-2:4-2:5 (2x1 + 2x2 + X3 + 3x5)! ' 



(2.5) 
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The number of ways to choose the leaves of T5 is 



n/ 2x1 + 2a;2 + X3 + 8x5 - 3fc\ _ (2a;i + 2x2 + + 8x5)! 
^^^j V 3 j " (3!)-5 (2X1+2X2 + X3)!' 

The number of ways to choose the leaves for the completion of T3 is 

""fj^ (2x1 + 2x2 + X3 - A;\ _ (2xi + 2x2 + X3)! ,^ 
1 )~ (2x1 +2x,)! • ^ 



fc=0 



For completing r2 and Ti , the number of ways to choose the leaves are 
given consecutively 



X2 — 1 



k=0 

and 



yr / 2X1 + 2X2 " 2k\ ^ (2xi+2x2)! 



(2!)^2 (2xi)! 



A;=0 ^ 

So far we have determined the graph corresponding to 5*2 but not chosen 
the pairs of points corresponding to its edges. The number of choices for 
these points is 

_ -^^3n/4:-X2-2X3-3X4 _ 2"j'^/^-X2-X3-X4-Xi _ ^'^n/2-2xi-X2-X3-2x4-X5 

X (d-4)^2+a;B (ci_5)^5_ (2.10) 

Having 5*1 and 52, we observe that in Ti there are 2xi vertices of degree 

3, while in T2 there are X2 vertices of degree 1 and X2 vertices of degree 5. 
In T;^ we have 2.T3 vertices of degree 1 and X3 vertices of degree 4, while in 
T4 we have 8x4 vertices of degree 1 and X4 vertices of degree 8 and in T5 we 
have X5 vertices of degree 6. For the isolated 3-stars in Si and ^2 we have 
2(ra/4 — X2 — X3 — X4 — X5) = ra/2 — 2x2 — 2x3 — 2x4 — 2x5 vertices of degree 

4. The remaining vertices (there are n/2 — X3 — 2x4 + X5 vertices) are of 
degree 2. Thus the number of free points in V is 

2{d - 8)xi + {d- l)x2 + (d - 5)x2 + 2{d - l)x3 + (d - 4)x3 + 3{d - l)x4 + {d- 8)x4 
+ {d- 6).T5 + {d- 4)(n/2 - 2x2 - 2x3 - 2x4 - 2x5) + {d- 2)(n/2 - X3 - 2x4 + X5) 
= (d - 3)n + 2x1 + 2x2 + 4x3 + 6x4. 

Therefore the number of ways to complete the pairing P is 

N{{d - 8)n + 2x1 + 2x2 + 4x3 + 6x4). (2.11) 
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Multiplying equations ( |2.3[ ) - ( |2.11 ) by the number of ways to choose Si 
as in (2.1), then dividing by N{nd), we have 

Ey;(y; - 1) = (^— j - 1) ((d - 2)(d - 3))^/^ 



X 
X 



y-r ^ (3n/4 — j;i — X2 — 2x3 — 3x4)!^ 



(n/4 - xi - X2 - X3 - X4 - xs)!^ (n/2 - X3 - 2x4 + X5)! 
32^1+2x2+2x3+^4 _ ^ 2x1 + 2x2 + 4x3 + 6x4)! 



2xi+x'2+x3+2x4 Xi!x2!x3!x4!x5! ((^^) n + xi +X2 + 2x3 + 3x4)! 

{d - 4)^2+X5 _ 5)X5 

(^(f _ ]^ja;2+2x3+3x4 _ 2^a;2+X3+a;4+a:5 (^^^ _ 3^2xi+X2+a::3+2x4+X5 

where 

Rl = {(xi,X2,X3,X4,X5) I Xl,X2,X3,X4,X5 > 0, Xl +X2 +X3 +X4 + X5 < n/4}. 

Set 

Xl X2 X3 X4 X5 

p=— , g = — , ?- = — , s = — , t=—. 
n n n n n 

Then Stirling's formula gives 

X ^ t)(F(p, (7, r, s, t)r (2.12) 

R2 



where 



F{p,q,r,s,t) 



R2 = {{p, q,r, s,t) I p,q,r, s,t> 0,p + q + r + s + t< 1/4} 

/(3/4-p-g-2r-3s)2 / ((^) +p + g + 2r + 3s) 



/(1/4 -p-q-r-s-t)'^ /(1/2 - r - 2s + f (p) f (q) f {r) f (s) f (t) 

^2p+2q+2r+s 2(i-3+p+g+3r+4s _ ^jg+t _ 



_ l)g+2r+3s (^^l _ 2)q+r+s+t (^^ _ ^-^2p+q+r+2s+t 

with /(x) = x^' and 

, , 2(3/4 -p-g-2r-3s)2 ^ 

a{p,q,r,s,t) - 



(1/4 - p - q-r - s-ty (1/2 - r - 2s + t) pgrst 

Since by convention /(O) := 1, it can be seen that F is continuous in R2. 
Next we determine the main contribution of the sum which comes from the 
maximum of F in R2. The following three lemmas prove that the maximum 
of F is attained at 



/ 9 9(d-3)(d-4) 9(d-3) 3 {d - 3){d - A){d - 5) 

V16d' 16d(d- l)(d- 2)' 8d{d-l){d-2)' 8d{d - l){d - 2)' 16d{d - l){d - 2) 
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with 



(2d)i-'^/2 (2(d-3/2))'^-3/2 
((d-l)(d-3))V2 ■ 



Lemma 2.1 [Ig Lei F and R2 be as in (|2J2j). Then for d>4, 
9(d-3)(d-4) 9{d-3) 3 



(d-3)((i-4)(d-5) 



\l6d' 16d{d - l){d - 2)' 8d{d - l){d - 2)' 8d{d - l){d - 2)' 16d{d - l){d - 2) 
is the local maximum point of F in the interior of R2 with 

(2(i)i-'^/2 (2(d-3/2))'^-3/2 



Pi^-^d—max) — 

Moreover for 4 < d < 10, x^- 
interior of R2. 



((d-l)(d-3))V2 • 
X is the global maximum point of F in the 



Proof. First we look for all critical points of F in the interior of i?2- We 
set the partial derivations of InF with respect to p, q, r, s and t, equal to 0, 
resulting in five equations 

= 9(1/4- p - q-r- s- tf{d - 3 + 2p + 2g + 4r + 6s) 

-p{d-2,f{'i/A-p-q-2r-'isf, (2.13) 
= 9(1/4 - p - g - r - s - t)2(d - 3 + 2p + 2g + 4r + 6s)(d - 4) 

-q{d -l){d- 2){d - 3)(3/4 -p-q-2r- 2,sf , (2.14) 
= 18(1/4 -p-q-r-s- tf{d - ?> + 2p + 2q + Ar + 6s)^(l/2 -r-2s + t) 

-r{d - lf{d -2){d- 3)(3/4 -p-q-2r- 3s)^ (2.15) 
= 6(1/4 - p- g- r- s- tf{d - 3 + 2p + 2g + 4r + 6s)^(l/2 -r-2s + tf 

-s{d - lf{d -2){d- 3)2(3/4 -p-q-2r- 3sf, (2.16) 
= {l/4-p-q-r-s- tf{d - 4)(c/ - 5) - t{d - 2){d - 'i){l/2 - r - 2s + t). 

(2.17) 

After substituting ( |2.13[ ) to ( |2.14| ) and ( |2.15[ ) to ( |2.16[ ) we have 

{d-l){d-2)q 
(c/-3)(d-4) ' 

and 

1 

t= ■ 



(2.18) 



r(.-3+^(^^:illi^ + 2, + 4. + 6. 



(d-3)(d-4) 



2r - 3s 



(1/2 - r - 2sK (d - 3 + 'l;^^!!^ + 2, + 4. + 6. 



(2,19) 
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After substituting and to ^J^, and ^7l7\i we have 

three homogenous equations 

Pi(g,r,s)_^ P2{q,r,s) ^ ^ P3{q,r,s) ^ ^ 



Qi(g,r,s) 



Q2iq,r,s) 



Q3{q,r,s) 



where 
Qi(g,r,s) 

Q2(g,r-,s) 
Q3(g,?^,s) 



256(d - 3f{d - Af {{Ad^ - 20d + 28)q + (4^^ - 28d + 48)r + {6d^ - 42d + 72)s 
+ - lOd^ + 33(i - 36) r^, 
8(d-3)((i-4)V gi(g,r,s), 

512 (d-3)4 (d-4)4 r2 



and degrees of q in Pi{q, r, s) , P2{q, r, s) and P3(g, r, s) are 2, while degrees 
of r and s are 5. It is sufficient to look at their nominator parts, 



Pi{q,r,s) 


= 0, 


(2.20) 


P2{q,r,s) 


= 0, 


(2.21) 


P3{q,r,s) 


= 0. 


(2.22) 



After taking the resultant of and ( [2^ and of ([23o]) and ([232]) with 

respect to q, we have two homogenous equations, 

Ui{r,s)Ti{r,s) = 0, (2.23) 
U2{r,s)T2ir,s) = 0, (2.24) 

where 

Ui{r,s) = -331776rS2(d-2)2(2c^-3)2(fi2_5d + 7)2(d-l)4(d-4)^(d-3)^ 
(d-3)4 



f/2(r,s) 



9r2 



f/i(r,s) 



and degrees of r and s in Ti(r, s) and T2(r, s) are 6. By taking the resultant 
of Ti(r, s) and r2('^, s) with respect to r, we have 

= -4497760410984972288si^(d - l)^{d - 3)^^(4^ _ i)4 

{8d^s - 2Asd^ + IQsd - 3) V{s). (2.25) 

where V{s) is a polynomial of degree 15 (see [13j for details). It is easy to 
show that one of the feasible solution for (2.25) is 

* = 3 
^ ~ 8d{d-l){d-2)' 

By substituting ( |2.26D to ( |2.23[ ) and ( |2.24| ) we have 

T^(r,s*) A^i(r,s*) = 0, 
VF(r,s*) A^2(r,s*) = 0. 



(2.26) 



(2.27) 
(2.28) 
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where 



W{r, s*) = 8rd^ - 24rd^ + 16rd + 27 - Qd 



and degrees of r in Ni(r, s*) and N2{r, s*) are 6. From W{r, s*) = we have 
a feasible solution for r 



9{d-3) 
M{d-l){d-2) 



(2.29) 



By substituting (2.29) to (2.20) - (2.22) we have three new equations 



Z[q,r\s*) Si{q,r*,s*) 
Z{q,r*,s*) S2{q,r\s*) 
Z{q,r*,s*) S3{q,r\s*) 



0, 
0, 
0. 



(2.30) 
(2.31) 
(2.32) 



where 



Z{q, r*,s*) = 16qd^ - 9d'^ - iSqd"^ + 32qd + 63d - 108 



and degrees of q in Si{q,r* , s*), S2{q,r* , s*) and £'3(5, r*, s*) are 2. From 
Z{q, r* , s*) = we have a feasible solution for q 



, ^ 9{d-3)id-4) 
^ 16d(d- l)(d- 2)' 

By substituting ( |2.33D to ( |2.18D and ( |2.19| ) we have 



(2.33) 



P 
t* 
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16d' 

{d-3){d-A){d-5) 
16did - l){d - 2) 



Then we obtain 



/ 9 9((i-3)(d-4) 9(d-3) 



{d-3){d 



Vl6d' 16d{d - l)((i - 2) ' 8(i(d - l)(d - 2) ' 8d{d -l){d-2)' I6d{d 



-4)(d-5) 
l)((i-2) 



Note: See [13J for details of Pi (g, r, s), P2{q,r,s), P3{q,r,s), Ti{r,s), T2{r,s), Ni\ 
N2{r,s*), Si{q,r*,s*), S2{q,r*,s*), S3{q,r*,s*) and V{s). 

For d > 4, the Hessian matrix of F at :x.d-max is negative definite. Then 
Xd-max is the local maximum point of F in the interior R2 with 

, (2d)i-'^/2 (2(d - 3/2))'^-3/2 



r, s 



P(xrf_ 



((d-l)(d-3))V2 



To prove that x^-max is the global maximum point in the interior of R2 
for 4 < d < 10, we use the following procedure, 



(1) For each d, determine the feasible solutions of (2.25). If there are k 



feasible solutions then denote them by s^i, . . . , Sdk, 
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(2) For s = Sdi, 1 < « < ^, determine the feasible solutions of (2.23) 



(2.24). If there are /j feasible solutions for each i, then denote them 
by ran, . . ■,rdm, 

(3) For s = Sdi and r = r^ij, 1 < i < k, 1 < j < li, determine the feasible 
solutions of (2.20) - (2.22). If there are niij feasible solutions for each 

1 Qdijuiij 1 



i and j, then denote them by q^iji 



(4) For s = Sdi, r = raij, and q = qdijf, 1 < i < k , 1 < j < k, 1 < f < niij 
determine Pdijf and tdijf in (2.18) and (2.19), 



(5) Define Xdij f = { pdij f ,qdij f,rdij, sgj, tdijf) as another feasible solution 
for system ( |2.13| ) - (2.17) in the interior of R2, 



(6) Determine F{xd- 



Fi^dijf)- 



For each d, 4 < d < 10, we have k = l,li = l^ruij = 1. It means that 



we have only one other feasible solution for (2.13) - (2.17) in the interior of 
i?2. Because F{xdijf) < F{xd-max) for each d, 4 < d < 10, then Xd-max is 
the global maximum point of F in the interior of i?2- CH 

To study the behavior of F on the boundary of R2, we generalize the 
approach used by Garmo in the proof of (|5J, Lemma 12). First let x = 
{xi,X2, ■ ■ ■ ,Xr) and Uj = n2,j, . . . , Ur-,j) for r > 2. The In function 

is defined on the set of non-negative real numbers with, by convention, 
X InO = 0. 

Lemma 2.2 Let R := R2 be a closed set in W and let 5R be the boundary 
of R. Assume that every point in 5R is the endpoint of an interval in R\6R. 
Let fi{x) = bi + UiX^ for i = 1, ■ ■ ■ ,m, where bi and Ui are constant, such 
that fi{x) > for all i and all x G R\5R. Define F to be a function on R 
such that 

m m 

F{x) = go{x) + ^aigi{x) = go{x) + ^ aifi{x) In fi{x). 



2 = 1 



i=l 



with Oi < for i < tuq < m. Suppose that for every x & R, the directional 
derivative of go at x in any direction is bounded. Let xq € 6R such that 
fi{xQ) = for at least one i < tuq and fi{xQ) > for all mo < i < m. Then 
Xq is not a local maximum of F on R. 

Proof. See [2] or [13] for details. 



Lemma 2.3 Let F and R := R2 be as in (2.12). Then the maximum 
of F does not occur in 6R. 
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Proof. We define x = {p, q, r, s, t) and ^ = {vi,V2, ws, ^4, fs)- Following the 
notation in Lemma 12.21 we write 



InF(x) = 5o(x) + ^ai/i(x)ln/i(x) 



where 

5o(x) = (d - 3 + 2p + 2g + 4r + 6s) In (d - 3 + 2p + 2g + 4r + 6s) + {2p + 2q + 2r + s) In (3) 
+ {q + t) In {d - A) + {t)ln{d - 5) - {p + q + r + 2s) In (2) -{q + 2r + 3s) In {d - 1) 
- (g + r + s + t) In (d - 2) - (2p + g + r + 2s + t) In (d - 3), 



ai = 02 = 03 

64 = &5 = 0, ftg 



04 = 05 = og = os = ~1) 07 = ~2, Og = 2, 61 
, 67 = 1/4, &8 = 1/2, bg = 3/4 and 



cZ-3 



Ul = 


(1,0,0,0,0), 


U2 = 


(0,1,0,0,0), 


U3 = 


(0,0,1,0,0), 


U4 = 


(0,0,0,1,0), 


U5 = 


(0,0,0,0,1), 


U6 = 


(1,1,2,3,0), 


U7 = 


(-1,-1,-1,-1,- 


Us = 


(0,0,-1,-2,1), 


Ug = 


(-1,-1,-2,-3,0) 



1), 



For gQ we have 
d 



dc 



50 (x + cv) 



c=0 



{2vi + 2t>2 + 4f3 + 6f4) In (d - 3 + 2p + 2(7 + 4r + 6s) 

+ {2vi + 2v2 + 4w3 + 6V4) + {2vi + 2^2 + 2v^ + v^) In (3) 
+ {v2 + W5) In (d - 4) + (fs) In {d - 5) 

- (f 1 + t;2 + W3 + 2^4) In (2) - {V2 + 2^3 + 3174) In (d - 1) 

- {v2 + f3 + "^4 + v^) In {d-2) - {2vi + ^2 + ^^3 + 2f4 + W5) In {d - 3) 



which is bounded for all x € i?. Having og > 0, then from Lemma 2.2 
only need to consider the solution of the following system 



we 



/^(x) 

/9(X) 



> 0, i 
= 



1,... 
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which is equivalent to 



p > 0, 

q > 0, 

r > 0, 

s > 0, 

t > 0, 

d-3 

+ p + q + 2r + 3s > 0, 



2 

1/4- p - q - r - s - t > 0, 

1/2-r -2s + t > 0, 

3/4-p-q-2r-5s = 0. 

It is easy to show that the only solution for the system is ci = (0, 0, 0, 1/4, 0). 
Consequently for xq E Ji?\{ci}, InF and xq satisfies the hypotheses of 
Lemma 2.2 Hence F does not have any maximum on (5i?\{ci}. Moreover 
for d> 4, 

3^/^ (2(d-3/2))'^/^-3/4 

~ ((d- l)3(d-2)((i- 3)2)1/4 

is strictly less than F (xd-max) ■ 



In the following lemma we will show that the sum in (2.12) can be 
approximated within a small region around the maximum. 

Lemma 2.4 [13\ Let B = B{X(i_max-,^) be a ball centered at 

_ /_9_ 9{d-3){d-4) 9{d-3) 3 {d - 3)id - 4){d - 5) 

Xd-max - yj^, 16d(d- l)(d-2)' 8d{d-l){d-2)' 8d{d - l){d - 2)' 16d{d - l){d - 2) 



and diameter 6 := n ^1"^ . Then with F and R := R2 as in (2.12), we have 
Y,a{x)F''{x) ~ ^a(a;)F"(a;). 



R 

Proof. Write 



^ a(x)F"(x) = a(x)F"(x) + a(x)F"(x). 

R B R\B 



It will be shown that 



^ q(x)F"(x) = o{a{^d-raax)F'^{yid-max))- 
R\B 
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For X G i?, the Taylor expansion of F at ^d-max is 

(5 5 
i=l j=i 

where 

9 9(d-3)(d-4) 9(d-3) 
■5l=P-T7ri, 52 = g - -.^ „ , TV73 53 = ?' 



16d' ' ^ 16d(d- l)(d-2)' ^ 8d(d- l)((i-2)^ 

3 _ (d-3)(rf-4)(d- 5) 

~ ^~ 8d(d- l)(d-2)' 16d(d- l)(d-2) 

and Cij are coming from Hessian matrix, H = (aij) G Ms^s, of -F in the 



proof of Lemma 2.1 



2 (^ij ) if ^ — j 

Uij, ifi/i 

For X* G where is the boundary of i?, we note that 

0{e^-''') = o(l). 

Then 

a(x*)F"(x*) ~ a{xd~max)F''{^d-max)o{l) = o{a{^d-max)F'' {^d-max)) ■ 

Since F attains its maximum uniquely at Xd-max, then for x G R\B 



a(x)F'^(x) = max a(x*)F''(x* 



Thus we have 



^ a(x)F"(x) = o{a{xd-max)F"'{xd-max))- □ 

Now we determine Y2b ct(x)i<'"(x). Since the summation concentrates 
near the maximum, each term a{x.) can be taken as a{xd-max) with 

8192^6 - l)3/2(d - 3/2)((i - 2)2 

O^y^d—max ) — 



243 ((i-3)5/2(d-4)((i- 5)1/2 



Referring to the Taylor expansion of F as in the proof of Lemma 2.4 
have 



we 



5 5 



^a(x)F"(x) a{xd-max)F"'{^d-max) ^^^P ~^ XI ^'^ij^i^j 
B B \ i=l j=i 
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The summation is a Riemann integral for the five-tuple integral 

5/2 



n 



/n fix i-n i-n i-n I 5 5 \ 

j j j j 6xp — n Cijtitj I dti dt2 dt^ dt^ dt^ 

-n J —n J —n J —n J —n \ »■ i „ / 



where 



n 



^2 



9{d-3){d-4) 
^ 16d{d-l){d-2) 



n 



9{d - 3) 



8d{d-l)id-2] 



n 



n 



8d{d-l){d-2] 



n 



t5 



^ (d-3)(d-4)(d-5) - 
l6d{d-l){d-2) ' 



n ^/n 
As n — )• oo, the range of the integration can be extended to ±00 without 
altering the main asymptotic term. Thus it is asymptotic to 

5 5 
exp I - I X] Yl 



n 



5/2 



00 roo foo foo foo 



00 J ~00 J ~00 J — 00 J — CO 



dti dt2 dts dti dt^. 



1=1 j=i 



The evaluation of the integral results in 

162^6 7r5/2(d - 3/2){d - 3f/\d -A){d- 5)^^ 
d^{d -l){d- 2)2(4^3 - 13(^2 + 36(i _ 36)V2 

Then (2.12) becomes 

■(cZ-2)(d-3)V/^' 



Ey;(y;-i)~n5/2x 



8(n7r)5/2 



V2 (2(i) 1-^/2 



3! 



'(2d)i-^/2 (2(d-3/2))"'-3/2\" 8192^/6d3(d-l)3/2((^_3/2)(d-2) 



((d_l)(d-3))i/2 



243 (d-3)5/2(d_4)(ci_ 5)1/2 



162^/6 ^5/2(tf _ 3/2)(d - 3)3/2(d - 4)(d - 5)^/2 
d^{d -l){d- 2)2(4(i3 - 13d^ + 36d - 36)1/2 



Kd- l)i/2(rf_ 3/2)2 



(d - 3)(4d3 - 13^2 + 36d - 36)i/2 



3/2)^-3/2 



(d- l)(d-2)V/2 



3! 



(2.34) 



As BY^ ^ 00 for d > 4 we have Eyj"(yj' - 1) ~ BY^"^. Thus from 



(2.34) and (2.1) we have 



EYf ~ 



2(d- l)i/2(d-3/2)^ 



^ (Ey,*)2. (2.35) 



(d - 3)(4d3 - I3d2 + 3Qd - 36)1/ 



Since VarY^ = BY^^-CEY^)^, the above equation gives the required result. 
□ 
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3 Expectation conditioned on short cycle distri- 
bution 



Lemma 3.1 \13\ Let n = mod 4 and be the number of cycles of length 
k in G{P) for P G Vn,d- Then for any finite sequences ji, . . . , jm of non- 
negative integers 



k=i 



where 



{d-1 



2k ' 

/ -3{d - 2) + V-15d2 + 2Ad\ ^ ( -3{d - 2) - V-I5d^ + 2M\ ^ 
^ ~ y 4((i- l)((i-3/2) j ^ y 4(d-l)(d-3/2) J ' 

Proof. To prove the lemma we first establish 

HY*Xk) I I ( -3,{d-2) + ^-lM^ + 2Ad V 

Ey; ""2^1^ \ 4(d-3/2) ] 

/ 3(rf_2)-V-l5Q;2 + 24(i y 
+ 4(d-3/2) ] ' ■ ^■^■^^ 

The number of ways to choose a cycle of length k in the pairing (with a 
distinguished point in a pair) is 

{d{d-l))K (3.2) 



(n- A;)! 

This induces an orientation and also a distinguished edge called a root 
edge in the cycle. 

Let C denote the set of pairs that corresponds to an oriented and rooted 
fc-cycle. Define S to be the set of pairs corresponding to a 3-star factor. 
Fix C and suppose C r\S consists of sq 0-stars (By 0-star we mean isolated 
vertices) lying at the centers of stars in the 3-star factor, si 1-stars, S2 2- 
stars and s^ = k — sq — 2si — 3s2 0-stars lying at the leaves of the 3-star 
factor. 

The edges of C can be classified into three types. The first type is the 
edges not lying in the 3-star factor, we denote this as 0. The second is the 
edges of 1-stars and the first edges of 2-stars and the last is the second edges 
of 2-stars. We denote them as 1 and 2 respectively. If we walk along C from 
the root edge, we obtain a sequence 5*0 E {0, i, 2}'^. 
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Fix C and Sq. The number of ways to choose the centers of the remaining 
(n/4 — Sq — s\ — S2) 3-stars, together with the points used, is 



Vn/4 - So - si - 52/ 



(3.3) 



The number of ways to choose the points in the centers of (sq + si + S2) 
3-stars is 

(((i-2)((i-3)((i-4))""((d-2)(d-3))^i(t^-2)'' = (d-2)^»+^i+^2(^_3)so+si(^_4^so_ 

(3.4) 

The number of leaves remaining for the 3-star factor is 

3(n/4 — Sq — si — S2) + 3so + 2si + S2 = 3n/4 — si — 2s2- 

The number of ways to select the leaves from the remaining (n/4— sq— si— S2) 
3-star is 

n/4-so_^i-s2-i ^^^^ _ ^ (3n/4-si-2g2)! 

l}^ V 3 y (3!)'^/4-*o-*i-*2(3so + 2si + s2)!' 

(3.5) 

The number of ways to choose the leaves from si 1-stars 

'y\ f^^o + 2si + S2 - 2k\ _ (3.S0 + 2.S1 + .S2)! 

l[\ 2 (2!)-(3so + .2)!- ^'-'^ 

The number of ways to choose the leaves from S2 S-stars 

'ri ( •^'^0 + S2-k\ _ (3so + g2)! „s 

The number of ways to choose the leaves from sq 0-stars 

n'°"Y3so-3A;\ (3so)! 

...oV 3 J-(3!F- ^'-'^ 

The number of ways to choose which vertex will be the center of s\ 

2*1. (3.9) 

The number of leaves of 5 lying in the cycle is A; — sq — 2si — 3s2- Note that 
every vertex in the cycle uses their two points. Thus the number of ways to 
choose the points that represent the leaves is 

(d-2)'=-*o-2*i-=^*2_ (310) 
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There are {n — k) — (n/4 — sq — si — S2) = (3n/4 — k + sq + si + S2) leaves 
outside the cycle. Note that there are d points in every vertex of the cycle. 
Thus the number of ways choose the points that represent the leaves is 



^3n/4:-k+so+si+S2 



(3.11) 



By multiplying (3.5) - (3.11) we have the number of ways to choose the 
leaves (including the points used) 

(3n/4 - Si - 2S2)! _ 2jfc-so-2si-3s2 ^3n/4:-k+so+si+S2 



(3!)"/4- 



S1-S2 



(3.12) 



The number of points for the leaves of 3-stars outside the cycles is {d - 
l)(3n/4 — /c + so + si + S2), while for the centers is (d — 3)(n/4 — sq — si — S2) 
The number of points for the centers of 3-stars inside the cycle is {d — 3)k- 
2so + 2s2- Thus the number of free points in the pairing is 



(d - l)(3n/4 -k + so + si + S2) + id- 3)(n/4 
= n(d - 3/2) -2k + 2si + 4s2. 



•So - si - ■S2) + {d- 3)k - 2so + 2s2 



Hence the number of ways to complete the pairing is 
N {n{d - 3/2) -2k + 2si + 4^2) . 



(3.13) 



Multiply ([3^ - by ( |3.12D - ( |3.13D , sum over all possible Sq, then 
multiply by (3.2). This results in the number of pairings containing a 3-star 
factor and an oriented and rooted cycle 



n — k 
n/4 - So - si 



S2 



So 

X (d-2)-+-+- (d-3)-+-(d-4)^o (d 



(d(d- l)((i-2))"/^-"o-"i-"2 

so — 2si—3s2 

(3.14) 



^ ^3n/4-k+so+sj+S2 ^ (^(^ _ 3/2) -2k + 2si + 4S2) . 

Dividing ( |3.14[ ) by the number of pairings with a 3-star factor, which is 



n! 



(n/4)! 



((i-l)((i-2)V/^ 



3! 



and then evaluating asymptotically, we obtain 



3{d-l){d-2) 
4{d - 3/2) 



E 

So 



(d-3)(d-4) 
3(d- l)(d- 2) 



so 



iV(d(n- 3/2)) 



8{d-3/2){d-3) 
3{d-l){d-2f 



32{d-3/2f 



%d-l){d-2Y 
(3.15) 

We follow an approach used in [7] to determine the summation. We can 
view 0, 1, 2 as three states in a Markov Chain, where the final state is 
equal to the initial state. We observe that 
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• 1 followed by means we pass a 1-star and this contributes a factor 

8(d-3/2)(d-3) 
3(d-l)(d-2)2 ' 

• 1 followed by 2 means we pass a 2-star and this contributes a factor 

32(d- 3/2)2 
9(d- l)(d-2)3' 

Thus for the transition matrix given by 



/ ^ ^ (d-3)(d-4) ^ 



3((i- 2) 
8(d- 3/2)((i- 3) 

3(d- l)(d-2)2 " 9(d-l)(d-2)3 
10 









32(d- 3/2)2 



we have 



TrlA'^) = 

So 



(d-3)((i-4) V° /8(d-3/2)(d-3) 



32(d- 3/2) 



2 \ S2 



^3((i-l)(d-2)y V l)(d-2)2 ; V9(d- l)(f^-2)3 

Since the eigenvalues of A are 

4(d-3/2) 



then 



rr(A'=) 



7l 
72 
73 



4(d-3/2) 
3{d-2) 



3{d-2) ' 

-3((i - 2) + V-15d2 + 24d 
3((i-l)((i-2) ' 

-3((i - 2) - V-15d2 + 24d 
3((i-l)(d-2) ' 



+ 



-3{d - 2) + ^-15^2 + 24d \ 
3((i- l)(d-2) ) 



+ 



-3((i - 2) - ^-15^2 + 24^' 
3((i- l)(d- 2) 



Then (3.15) becomes 



(^-1)^^ + 



-3(d - 2) + V-15d2 + 24(i' 
4(d-3/2) , 



+ 



-3((i - 2) - V-15(i2 + 24d' 
4(d-3/2) 



(3.16) 

After (3.16) is divided by 2k to remove the orientation and rooting of the 
cycle we obtain (3.1). □ 
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4 Proof of Theorem 



1.1 



First we prove the following theorem. 

Theorem 4.1 Let n = mod 4 and A < d < 10. Then for P G Vn,d, G{P) 
a.a.s has a 3-star factor. Moreover, 



1^ = JJ(1 + Skf^'e-^"^^ forn^oo 



EY* 

d- k=l 

where are independent Poisson variables with EZk = \k for k > 1 and 
{d-1)^ 



2k ' 

-3{d-2) + V-l5d^ + 2mY ( -3(d - 2) - V-lSd^ + 24d\ ^ 



y 4((i- l)((i-3/2) I ^\ 4(d-l)(d-3/2) j 

Proof. We will show that satisfies the conditions (A.l) - (A. 4) in (The- 
orem 4.1, [H]). Since is the number of short cycles of length A; in a 

pseudograph coming from Vn,d-, then (A.l) is fulfilled with = ^'^'2^ , by 
Bollobas' result on short cycles in Vn.d [3j. The condition (A. 2), (A. 3) and 



(A. 4) are fulfilled consecutively by Lemma 3.1, Theorem 2.1 and Theorem 
[221 □ 



Proof of Theorem 11.11 

Theorem |1.1| comes directly from Theorem 4.1 by Lemma |1.1[ From the 
argument in (Remark 9.25, [6j) we also obtain 

BYd /3(5(i2 - l2d + 6) 

EI7 ^ V 4(2d - 3)2 

EFj / -9{8d^ - 63d^ + 206^^ _ 322^^ + 216d - 36) 



(Ey,)2 ' ^""^V 4(2d-3)4((i- 
f 2(d- 1)1/2(^-3/2)2 



y{d- 3)(4d3 - I3(i2 + 36d - 36)1/2 J ' 

We should point it out again, that most part of the computation in 
determining the second moment of the number of 3-star factors in G S Vn,d 
is valid for general d> 4. The only part that is still hard to prove is showing 
that the desired maximum point is the global maximum. From what we have 
in the case 4 < d < 10, we conjecture that the asymptotic distribution of the 
number of 3-star factor in random d-regular graph have the same behaviour 
for d> 4. 



□ 
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Conjecture 4.1 Restrict n to mod 4 and d > 4. Then G G Qn,d a.a.s has 
a 3-star factor. Furthermore, letting Yd denote the number of 3-star factors 
inG e gn,d, 

y. °° 

—^r ^W= Y\il + 6kf''e'^''^'^ forn^oo 
Eld 

where are independent Poisson variables with EZ^ = \k for A; > 3 and 
_ {d-lf 

/ -3(d - 2) + V-15d2 + 24d\ ^ ( -3{d - 2) - V-I5d^ + 24d\ ^ 
~ \^ 4(d- l)((i-3/2) J ^ 4((i-l)(d-3/2) J 
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